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Abstract 

A combination of finite energy sum rule techniques and Chiral Perturbation Theory (xPT) is 
used in order to exploit recent ALEPH data on the non-strange r vector (V") and axial-vector 
(A) spectral functions with respect to an experimental determination of the quantity 
Lio- A constrained fit of E^y_ A inverse moments (I < 0) and positive spectral moments 
(I > 0) adjusts simultaneously Lio and the nonperturbative power terms of the Operator 
Product Expansion. We give explicit formulae for the first k = 0, 1 and I = —1, —2 strange 
and non-strange inverse moment chiral sum rules to one-loop order generalized xPT. Our 
final result reads L r w (M p ) = -(5.13 ± 0.19) x 10" 3 , where the error includes experimental 
and theoretical uncertainties. 
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1 Introduction 



The nonperturbative features of strong interactions make QCD a rich environment for 
theoretical investigations. At sufficiently high energies it is possible to parametrize 
the nonperturbative effects by vacuum condensates, following the rules of Wilson's 
Operator Product Expansion (OPE) The universal character of these condensates 
has been used in the derivation of the so-called QCD spectral sum rules allowing, 
in principle, their determination from experiment. A particular role is played by the 
condensates which are order parameters of the spontaneous breakdown of chiral sym- 
metry (SB%S). The latter vanish at all orders of perturbation theory and they control 
the high energy behavior of chiral correlation functions, such as the difference of vector 
and axial current two-point functions. On the other hand, at low energies, SB^S makes 
it possible to construct an effective theory of QCD, the Chiral Perturbation Theory 
(xPT) [[| f|, which uses the Goldstone bosons as fundamental fields and provides a 
systematic expansion of QCD correlation functions in powers of momenta and quark 
masses. All missing information is then parametrized by low-energy coupling constants, 
which can be determined phenomenologically in low-energy experiments involving pi- 
ons and kaons. The fundamental parameters describing chiral symmetry breaking, 
the running quark masses and the quark anti-quark condensates (qq) appear both in 
low-energy (%PT) and the high energy OPE expansion. For this reason it is useful to 
combine the two expansions in order to get a truly systematic approach to the chiral 
sum rules ||. In this paper the combined approach is illustrated through a determi- 
nation of the Lio constant of the chiral lagrangian, including high-energy corrections 
coming from the OPE. The connection between the two domains is provided by ex- 
perimental data on r hadronic spectral functions published recently by the ALEPH 
Collaboration || |7j. 

At the leading order of %PT, Lio is directly linked to the vector, vi, and axial- vector, 
ai, spin-one spectral functions (the subscripts refer to the spin J of the hadronic sys- 
tem) through the Das-Mathur-Okubo sum rule [Q 

j-j J ds- [vt{s) - a x {s)] ~ -4L 10 . (1) 
o 

As it stands the DMO sum rule (|l|) is subject to chiral corrections due to non-vanishing 
quark masses ||. On the other hand, the integral has to be cut at some finite energy 
sq < M%, since no experimental information on v\ — a x is available above M%. This 
truncation introduces an error which competes with the low-energy chiral corrections. 
Both types of corrections can be systematically included through i) the high-energy 
expansion in a s (so) and in inverse powers of so, and ii) the low-energy expansion in 
powers of quark masses and of their logarithms. 
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2 Spectral Moments 



Using unitarity and analyticity, the spectral functions are connected to the imaginary 
part of the two-point correlation functions, 



n^(g) = l Jd 4 xe^(o\T(u? J (x)u»(ono) 



= y + q , qV) n (i )[/(g2) + tft n (o^ (g2) j (2) 

of vector (U£ = = qj^qi) or axial-vector (L^- = = <?j7 M 75<?i) colour-singlet 
quark currents for time-like momentum- squared q 2 > 0. Lorentz decomposition is used 
to separate the correlation function into its J = 1 and J = parts. The correlation 
function (Q) is analytic everywhere in the complex s plane except on the positive real 
axis where singularities exist. Using the definitions adopted in Refs. || |7J together 
with Eq. (0), one identifies for non-strange quark currents 

Due to the conserved vector current, there is no J = contribution to the vector 
spectral function, while the only contribution to ao is assumed to be from the pion 



pole. It is connected via PCAC to the pion decay constant, a 0:7T (s) = 47r 2 /^ S(s 



m 



2 ■ 



According to the method proposed by Le Diberder and Pich |L0], it is possible to 
exploit the information from the explicit shape of the spectral functions by calculating 
so-called spectral moments, i.e., weighted integrals over the spectral functions. If 
W(s) is an analytic function, by Cauchy's theorem, the imaginary part of , A is 
proportional to the discontinuity across the positive real axis: 

SO 

I ds W(5)Imn^ v/A ( S ) = -j i fds W{8)U$ y/A ( 8 ) , (4) 

|s|=so 

where so is large enough for the OPE series to converge. The authors of [|HJ choose 
for W{s) the functions 

W^\s) = (l-^) +k (-1 , (5) 



so/ vso- 

with k and / positive integers. The factor (1 — s/s Q ) k suppresses the integrand at the 
crossing of the positive real axis where the validity of the OPE is questioned. Its coun- 
terpart (s/so) 1 projects on higher energies. These moments were successfully applied 
in order to constrain nonperturbative contributions to the r hadronic width, R T , a 
procedure which lead to precise determinations of a s (M 2 ) 0, |TT], [12 . 



The extension of the spectral moment analysis to negative integer values of I ( "in- 
verse moment sum rules", (IMSR) flI5| ) requires, due to the pole at s = 0, a modified 
contour of integration in the complex s plane, as shown in Fig. [I]. This is where ;\PT 
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comes into play: along the small circle placed at the production threshold, s t h = 4M 2 
for non-strange (ud) and s t h = (M n + M K ) 2 for strange (us) currents, we can use xPT 
predictions for the two-point correlators. Using the weight function (^j we adopt the 
following definition of the moments: 



as I s \ "" / s 



Ml 



2+k / \l 



1 + 2 



imn 



(o+i) 

V/A 



2^Imn^( S ) 



(6) 



where s m i n = for the positive moments[] (/ > 0) and s m i n = s t h, which is the continuum 
threshold, for the inverse moments. According to the relation (f|), Eq. (Bf) reads 



R 



(k,i) 

t,V/A 



67ri|V u d| 2 5 , EW 



c 



M T 2 ^ M T 2 J 



2+k 



Ml 



x 



"( 1 + 2 ^)4«»( s )-2^n«»', 



(7) 



where C = C\ + Ci for the inverse moments and C = C\ for the positive moments (see 
Fig- B). 

1 This is due to the pion pole which is at zero mass in the chiral limit. 
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Due to the cut of the integral (||) at M 2 , nonperturbative physics parametrized by 
the short- distance OPE for scalar operators |], [| must be considered: 



n 



V/A 



»= E 



E 4m(«.M 2 )<Owa(a« 2 )> • 



0=0,2,4,. 



(8) 



dimO=D 



The parameter /i separates the long-distance nonperturbative effects, absorbed into 
the vacuum expectation elements (CV/a(/^ 2 )), from the short- distance effects which 
are included in the Wilson coefficients Cv/a(s, fi 2 ) M. We will assume the convergence 
of the OPE series at the r mass. This is justified in the light of the success of the 
analysis performed in Ref. |7j (see Ref. |15|] for details). Using the formulae of Refs. |14 



and [|T6j for the nonperturbative power expansion of the correlators, one obtains for 
the (V — A) case 



n (0+1) (s) 

LL ud,V-A\ b ) 



a.Js m'(s) 



7l z 



59 



rh{uu + dd) 16 m 4 (s) 



77T 2 S 2 



87r 2 a s (/x 2 
2tt 2 



— - -L(s) ) a s (n 



<<W)> 



4 ' 



hSv-a(-«) 



7T^ 



2al 



5 + 



21373 75 , , 

+ T77C(3) a s (s) 



rh{uu + dd) 



2448 34 
53 



m 2 (s] 



(9) 



2 m 2 (/x 2 ) 



7.U2 12 »-'« 



(10) 



with a s (s) 



a s (s) /tt, L(s) = log(s//i 2 ) and the dimension D = 6 operators 



C 1 

^6 



o 2 



. u-f^ 5 T a dd^-f 5 T a u - u ltl T a dd^T a u 
= uj^dd^u - U7 M 7 5 dd7 M 7 5 u, 



ir 



where the S'C/(3) generators T a are normalized so that tr(T a T b ) = 5 ab /2. We use the 
average mass m = (m u + rrid)/2 in the above equations, i.e., we assume SU(2) sym- 
metry. The constant C(fi 2 ) depends on the renormalization procedure^] and should not 
affect physical observables. The dimension D = contribution is of pure perturbative 
origin and is degenerate in all-orders of perturbation theory for vector and axial- vector 
currents. Dimension D = 2 mass terms are calculated perturbatively to order a 2 
which suffices for the light u, d, quarks. The coefficient functions of the dimension 
D = 4 operators for vector and axial-vector currents have been calculated to sublead- 
ing order in Refs. |T7|, |TB| . Their vacuum expectation values are expressed in terms of 
the scale invariant gluon and quark condensates. Since the Wilson coefficients of the 
gluon condensate are symmetric for vector and axial-vector currents, they vanish in 



2 We will assume a renormalization scheme that preserves chiral symmetry, so that C is the same 
for the vector and axial correlators. 
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the difference. The expectation values of the dimension D = 6 operators (|TTD obey the 
inequalities (0\(ii 2 )) > and (Og(/i 2 )) < 0, which can be derived from first principles. 
The corresponding coefficient functions were calculated by the authors of Refs. |19| in 



the chiral limit for which the J = contribution vanishes. For the dimension D = 8 
operators no such calculations are available in the literature, and we will assume that 
there is no logarithmic s dependence in leading order a s . Again, the J = contribution 
vanishes in the chiral limit. 

As constraints on the nonperturbative phenomenological operators introduced in 
Eqs. (|9]) and ( |T0D from theory alone are scarce, we will benefit from the information 
provided by the (V — A) spectral moments in order to determine the magnitude of 
the OPE power terms at M 2 . We therefore perform a combined fit of the IMSR (i.e., 
I = — 1) which determines L 10 , and the / > moments which adjust the nonperturba- 
tive contributions. 



3 Chiral Perturbation Theory 

The non-strange correlators (Q) have been calculated at one- loop level [|], 20 1 and, most 



recently, at two- loop level [21, |22[] in Standard %PT. In this paper we stick to the 0(p 4 ) 
one-loop order for the following two reasons: i) the high energy corrections are often 
more important than the 0(p 6 ) chiral corrections (whose precise estimate has not yet 
been fully completed |§) and ii) it is important to proceed in the combined analysis 
order by order in quark masses. On the other hand, we use the generalized version of 
%PT (GxPT) p3| , which allows us to investigate the sensitivity of the analysis to the 
variation of the quark condensate and of the quark mass ratio r = m s /rh. The standard 
XPT assumes |2j 2m(qq) ~ -F 2 M 2 and r ~ 2M|-/M 2 - 1 » 25.9, whereas G^PT 



admits lower values of these two quantities p5j, P3| . It is interesting to investigate 
whether the ALEPH spectral function data are precise enough to have any impact on 
the on-going debate about the size of (qq). Anyhow, the alterations of the standard 
0(p 4 ) results for non-strange correlators (|2|) introduced by G^PT are marginal. They 
merely concern the symmetry breaking J = component of the spectral functions and 
most of them are actually absorbed into the renormalization of F n (Fk). 

In order to make our analysis as independent of a particular truncation of the %PT 
series as possible, we proceed in two steps. First, one defines a phenomenological 
quantity called Lf^ via the contribution of the small circle C2 (see Fig. [I]) to the 
integral (^) of the chiral combination V — A for I = — 1. L|q is then determined in the 
combined fit of the IMSR and / > moments. The result of this fit is independent 
of the ;\RT renormalization scale /x x pt- The latter is used in the next step in order 
to relate Lf^ to the quark-mass independent, scale dependent constant Z^ (/i x PT) and 
finally to other observables (from ir —>■ ewy data, (r 2 ) 7r ). 
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3.1 Non-strange IMSR's 

For the Cabibbo-allowed channel we obtain: 

V$$(s) = 4M^( S ) + 8M^(,)-4(^ + 2iJ[), (12) 

nSvto = o. ( 13 ) 

= -^^-4(2^[-Ly , (14) 

= -^^ + 8m 2 (H 2j2 -2B 3 ), (15) 

while, for the strange channel the correlators read 

iCvW = GM r Kv (s) + 6M r Kw (s) - A(L[o + m) , (16) 

sU^ v (s) = Q(L Kv (s)+L K7T (s)) + 2m 2 (r-l) 2 (2B 3 + H 2 , 2 ) , (17) 

n£!?(*) = ~ 4(2i?[ ~ Lg ' (18) 

9 F 2 A/f 2 

Sn£) > } = + 2- 2 (r + l) 2 (^, 2 -25 3 ), (19) 



The functions Mp p/ (s) and Lpp>(s) are loop integrals, defined, e.g., in Ref. pOfl . The 
superscript r refers to renormalized quantities, which depend on the scale // x pt- The 
whole expressions are /x x pt independent. H 22 and £? 3 are found to be finite, in agree- 
ment with PB[ , and do not need renormalization. £f[ and H 2 2 are coefficients of contact 
terms of the sources. They are counterterms needed to renormalize the ultraviolet di- 
vergences of the Green functions and do not appear in physical observables. Our aim 
is to determine L\q: therefore we will consider the difference between the vector and 
the axial-vector correlators for which the constant Hi disappears. Correspondingly, as 
already pointed out, we will not need the perturbative expressions which are identical 
for vector and axial-vector cases. As for the constant H 22 which multiplies the term 

of the £(2,2) chiral lagrangianQ , it always appears in the same combination with C*(/i 2 ), 
in such a way that the ambiguities cancel out. We thus define a -£^2,2, in which the 
constant C{M 2 ) is absorbed. What is new at this order with respect to S%PT is the 
appearance of the constant B 3 which multiplies the term 

(tfD^D^x + h.c.) 



£(n,m) collects terms in the chiral lagrangian with n covariant derivatives and m powers of quark 
masses. In the same notation the Hi constant introduced by Gasser and Leutwyler |po|] would become 
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of the £(2,2) lagrangian. As can be seen from its form it is difficult to find a process 
in which B 3 would contribute directly. It will contribute to off-shell vertices involving 
Goldstone bosons. 

The non-strange IMSR's corresponding to I = — 1 and k — 0, 1 read: 



1 



|VU 2 & 



-R 



ud\ ^EW 



(0,-1) 
t,V-A 



p2 71^2 

96 I r 2 Lf ' + 24 I r 2 ' • 



T 

96 

+ J^Tr 2 m(dd + MM) 



M 4 

1 23 



12 



4896 + 34« 3 >H M '> 



17 



1 + + 



+ 



192vr 4 

4007T 4 

3M5 



a s (// 



/103 
\~24~ 
12 



- -L(M 2 )) a s (fi 



576 
7M 4 

(<W)> 



m 4 (M 2 ; 



a s (M 2 



a ' (/i2) ( 1 + i L(M ' } ) (0eV)) 



29 
24 



(20) 



1 

\Vud\ 2 SEW 



-R 



(i-i) 



- 967^ 



2 reff 
10 



24?r 2 — — 3 



F 2 M 2 /• 2 .U r 



i i " 7T _|_ T " " ; I 



144 

+ ^m 2 (M 2 ; 



+ 



M 2 
144 



M 4 .Uf 
71 /vr 2 39461 75 



+ 



M 4 
480tt 4 

M r 6 
472tt 4 



7r rh{dd + mm) 



a s (M 2 ) 24 
2 



1 12 IS*)Ci 31 



1 + -a s (M 2 ) 



43 



a 2 (M 2 



+ ^C(3) I a,(M 2 ) 

864 a4 2 , 
m (M T/ 



a s (/x 2 ) 



7M 4 
(0 6 V)> 



3M6 

where we have defined 



60 

59' 



W) <<W)> + 24tt 



a s (M 2 ; 



(21) 



8L 



cff 
10 



2s \ 



IFl 



F 2 

s - M 2 Ml 



Ml 



2 ) 

T / ' 



(22) 



which is proportional to the contribution of the small circle C 2 to the integral ([!]), with 
the pion pole subtracted. This quantity is a well defined observable, the ambiguity in 
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the two-point function being absorbed by the constant C. In the particular case of the 
one-loop G%PT calculation, its expansion reads: 

< - ^lofePT) + ^ (>o g g + l) + ^ log M + g ( 2fl3 - d w ) , (23) 

which is independent of /^ x pt- Unless stated otherwise all condensates, quark masses 
and xPT constants in the above expressions are evaluated at QCD renormalization scale 
Hqcd = M T , while the product of the light quark mass and the scalar quark operator, 
m(dd+uu) , is scale invariant. Taking the difference of Eq. fl20D and (|2lD and subtracting 
the contribution from the pion pole recovers the expression for R T y—A — R\- \/—a given 
in Due to the strong intrinsic correlations of 98% between the IMSR's defined 

above only one IMSR is used as input to the combined fit. We find it convenient to 
use the moment k — 1,1 — — 1 (Eq. (|2ID) because its experimental value is known with 
a 30% better precision which is due to the additional (1 — s/M%) suppression of the 
less accurate high energy tail of the (V — A) spectral function. 



3.2 Strange IMSR 



Analogous IMSR's with I = —2 would require the two-loop results for the correla- 
tors (g), because quark- mass independent terms from the £(e,o) lagrangian would give 
rise to new contributions which are not a priori small. However, if we consider the 
difference between the strange and non-strange (V + A) moments these terms cancel 
out because of the SU (3) symmetry. Of course there will be the two-loop corrections to 
the terms already present in the one-loop results, but these are subleading. Therefore, 
we can write a particular combination of inverse moments which does not contain any 

(1 —1) (2 —2) 

unknown low-energy constant. An example of such a combination is AR T V+A + R T V+A , 
so that the strange IMSR takes the form: 



S 



EW 



~ fi (4-^5 1 + R r,S 2 ) - \ v |2 i^ Rl r]v+A + R t,V+a) 
us * ud ' 



- (M^(0) + M^(0)) - Mj^(0) - 2ilC(0) 




+ 12 (M^(0) + M^(0)) - 8 (M r KK (0) + 2MIM) ~ ^ (L' Kt ,(0) + L^(0)) 
12 



12 /„ 



2 P 2 
K ^ n 



FX] 



s-7.91 



1) 



59 
15 



+ 



1927T 2 — j(m s ss — muu) 

My ' 

576 rh A 



7 M T 4 



a s (M T 2 



215 7 2 47 4 
j f T 

48 2 48 
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(24) 



where —7.91 is the value of the contribution, using the 7r° mass and the QCD K + 
mass (i.e., without electromagnetism) . We have neglected all logarithmic s dependence 
in the dimension D = 6 and D = 8 operators Oq AS ^ and Og AS \ The latter are expected 
to be suppressed because, in contrast to the non-strange C 6 and 0$, they vanish in the 
chiral limit. Data for the inclusive vector plus axial- vector strange spectral function 
from r decays are not available at present. Such data could provide information on the 
size of the quark condensate, which could represent up to 10% of the %PT contribution 
to the r.h.s. of Eq. 



4 Theoretical parameters and uncertainties 

(k I) 

When fitting the theoretical prediction of the R T y_ A moments to data, theoretical 
as well as experimental uncertainties and the correlations of these between the (k, I) 
moments must be considered. The masses of the light quarks are parametrized using 
the mass ratio r = m s /m of which the central value is set to the S^PT value of 26. 
A lower limit is found at r > r ^ m -,t = 2(M^/M vr ) — 1 ^ 6.1 (while rii m ; t ~ 8.2 when 
including higher orders p7[ ) which determines the range 

8 < r < oo . 

The average light quark mass is then obtained via rh = m s /r where we use for the 
strange quark mass m s (M T ) = 172 MeV/c 2 [28]. This parametrization makes it possible 



to use the theoretical correlation between rh and the quark condensate, which to leading 



order in quark masses is given by the generalized Gell-Mann-Oakes-Renner relation [25 



rn { uu + dd) ~ -FlMl K ' ^ ^ J , (25) 



(r — ri)(r + r\ + 2) 
r 2 — 1 

where r\ ~ 2(M^-/M 7r ) — 1. For the standard value r = 25.9, Eq. ( p5|) becomes the 
usual PCAC relation rh(uu + dd) = —F%M%. Corrections to Eq. ( [25] ) are expected to 
be small in the whole range of r so that we assume a relative uncertainty of 10%. We 
will comment in Section |6| on the sensitivity of the data with respect to the r ratio. 
Theoretical uncertainties are introduced from the strong coupling constant where, in 
order to be uncorrelated to the r data used in this analysis, we rely on the result from 
the global electroweak fit found recently to be p9l BUI 



a s (M 2 z ) = 0.1198 ±0.0031 . 

Uncertainties from the OPE separation scale /i are evaluated by varying [i from 
1.3 GeV to 2.3 GeV, while in the fit we choose fi = M T so that the logarithmic scale 
dependence of the dimension D = 6 terms vanishes after the contour integration. Ad- 
ditional small uncertainties stem from the pion decay constant, F n = (92.4 ±0.2) MeV, 
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taken from Ref. and the overall correction factor for electroweak radiation, Sew — 
1.0194, obtained in Ref. |j32| , with an estimated error of ASew = 0.0040 according to 
Ref. p|. 



An overview of the associated uncertainties in the theoretical prediction of the mo- 
ments is given in Table |I[ The moment errors from the a s uncertainty depend on the 
central input values of the nonperturbative operators. The numbers given in the fourth 
line of Table [L] correspond to the fit values, Eqs. (|30|)-(|31"|), which have been obtained 
in an iterative procedure. 



5 Spectral Functions from hadronic r decays 

The ALEPH Collaboration measured the inclusive invariant mass-squared spectra of 
vector and axial-vector hadronic r decays and provided the corresponding bin-to-bin 
covariance matrices || 0j. The mass distributions naturally contain the kinematic 
factor of Eq. (H) so that the measured spectral moments read 




1 

K 



(26) 



with the normalized invariant mass-squared spectra (l/Ny/A)(dNv/A/ds) of vector and 
axial- vector final states, the electronic branching ratio (using universality) |3T], |7], 
B e = (17.794 ±0.045)%, and the inclusive branching ratios 0, By = (31.58 ± 0.29)%, 
B A = (30.56 ± 0.30)%, as well as their difference, B V _ A = (1.02 ± 0.58)%. Due to 
anticorrelations between vector and axial-vector final states, especially for the KK7r 
modes where the vector and axial- vector parts are unknown, the error of the difference 
is larger than the quadratic sum of the errors on V and A. Figure |2] shows the (V — A) 
mass-squared distribution, which is the integrand of Eq. (|26D for zero moments, k = 
I = 0. With increasing masses it is dominated by the p (V), a± (A) and the p(1450), 
uj-k (V) resonance contributions which create the oscillating behaviour. Table [3] and § 
give the experimental values and uncertainties for the IMSR R T y_A and the k — 1, 
I = 0, ... ,3 moments as well as their correlations which are computed analytically from 
the contraction of the derivatives of the moments with the covariance matrices of the 
respective normalized invariant mass-squared spectra. 

Based on isospin invariance, the conserved vector current hypothesis (CVC) relates 
vector hadronic r spectral functions to isovector cross section measurements of the 
reaction e + e~ — > hadrons. There exist precise data on the low energy, time-like pion 
form factor-squared [^(s)) 2 measured by the NA7 Collaboration [34|. Using the CVC 
relation 

Vl,TT-no{s) = 




4ilf 2 \ 3/2 



^ =1 (*)| 2 , (27) 



10 
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> 

T3 
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Figure 2: Vector minus axial-vector (V — A) invariant mass-squared distribution mea- 
sured by ALEPH 0. 



one can include the additional data in order to improve the precision of the mo- 
ments (p6|), in particular for the IMSR in which the low-energy region is emphasized. 
Figure |] shows the vector spectral function from r data (three bins) together with the 
NA7 measurements for energy-squared s < 0.2 GeV 2 . In addition, we give the result 
when fitting both data sets using the parametrization 



FJs) 



\+ (r 2 ) n s + As 2 + Bs* 
o 



:28i 



for the pion form factor. Using analyticity, the pion charge radius- squared, (r 2 ) 7T 



(0.439 ± 0.008) fm , is taken from an analysis of space-like data 



We obtain the 



fit results A = -(7.5 ± 1.1) GeV -4 and B = (62.5 ± 6.4) GeV -4 with X 2 = 0.6 for 
5 degrees of freedom. The correlation between A and B is absorbed in the diagonal 
errors given, so that both quantities can be handled as being uncorrelated. Replacing 
for the above energy interval 4M 2 < s < 0.2 GeV 2 the pure r data by a combination of 
r and e + e~ data represented by the analytical expressions (|27|) and (^), we obtain the 
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results given in the third and fourth line of Table |T]. A small improvement in precision 
of 11% is observed for the IMSR. 

The spectral information is used to fit simultaneously the low-energy quantity Lfy 
and the nonperturbative phenomenological operators. For dimension D = 6 we will ne- 
glect the contribution of 0$, which is suppressed by ot 2 s and, furthermore, is suppressed 
relatively to 0\ in the large N c limit. Therefore we will simply keep Oq = 0\(M 2 ) and 
the 0$ operator of dimension D = 8. 



6 Results of the fit 

The fit minimizes the \ 2 °f the differences between measured and fitted quantities 
contracted with the inverse of the sum of the experimental and theoretical covariance 
matrices taken from Table The results of the fit are for Lf^: 



Lf f = -(6.36 ± 0.09 cxp ± 0.14 thco ± 0.07 fit ± 0.06 OPE ) x 1(T 3 



(29) 



and for the nonperturbative operators: 

(O e ) = (5.0 ± 0.5 cxp ± 0.4 thco ± 0.2 ftt ± Hope) X 10~ 4 GeV 6 , (30) 
(0 8 ) = (8.7 ± 1.0 cxp ± 0.1 t hco ± 0.6nt ± 2.1ope) x 10~ 3 GeV 8 , (31) 
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(k,l) ► 


(I, -I) 


(0,0) 


(1,0) 


(1,1) 


(1,2) 


(1,3) 


R [ T k f_ A (ALEPH) 


5. 16 


0.055 


0.038 


0.047 


-0.0164 


-0.0126 




0.09 


0.031 


0.017 


0.006 


0.0035 


0.0023 


R [ ry-A (ALEPH + NA7) 


5.13 


0.055 


0.037 


0.047 


-0.0164 


-0.0126 


A eXD DV^'V 


0.08 


n no i 

0.031 


0.017 


0.006 


0.0035 


0.0023 


A thco R ( r k p_ A (Ar) 


0.12 


0.003 


0.003 


0.001 


0.0003 


< 0.0001 




0.02 


0.009 


0.009 


0.002 


0.0029 


0.0001 


A theo R (h)_ A (ASew) 


0.02 


< 0.001 


< 0.001 


< 0.001 


0.0001 


< 0.0001 


Atheo^M^ ( A^ OPE ) 


< 0.01 


0.005 


0.005 


0.002 


0.0018 


< 0.0001 


Atheo^W^ (A (gg)) 


< 0.01 


< 0.001 


< 0.001 


< 0.001 


< 0.0001 


< 0.0001 


A thc 0jR W)_^ (Ai?7r) 


< 0.01 


< 0.001 


< 0.001 


< 0.001 


< 0.0001 


< 0.0001 


(Theory fitted) 


5.13 


0.061 


0.032 


0.053 


-0.0148 


-0.0098 



Table 1: Measured spectral Moments of vector (V) minus axial-vector (A) using r data 
only (ALEPH) and using r + e + e~ data (ALEPH + NA7). The quoted errors account 
for the total experimental uncertainties including statistical and systematic effects as 
well as the theoretical uncertainties according to Section [|. The last line gives the fitted 
theoretical moments using the parameters given in Eqs. (|2T?D-(PT|) . 



(M) 


(1,-1) (0,0) 


(1,0) 


(1,1) 


(1,2) 


(1,3) 


(i,-i) 


1 0.46 


0.61 


0.40 


0.26 


0.13 


(0,0) 


1 


0.89 


0.97 


0.84 


0.80 


(1,0) 




1 


0.88 


0.74 


0.45 


(1,1) 






1 


0.89 


0.78 


(1,2) 








1 


0.76 



Table 2: Sum of experimental and theoretical correlations between the moments R). v_a- 





Lio (O e ) (0 8 ) 


(Oe) 


1 -0.26 0.05 
1 0.14 



Table 3: Correlations between the fitted parameters fl2~9[)-(|31~|). 
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with a x 2 °f 2.5 for 3 degree of freedom. The errors are separated in experimental 
(first number) and theoretical (second number) parts, and a fit uncertainty (third 
number) is added. The latter is due to a well known bias when fitting quantities for 
which correlations are due to normalization uncertainties []36| (here the r branching 
ratios) leading systematically to lower values in terms of the normalization of the fitted 
parametrization. The errors quoted account for the differences between fully correlated 
and uncorrelated results. The authors of Ref. |7| observed a variation of the results on 
the nonperturbative operators depending on the weighting of the r spectral functions 
used in the actual fit. These variations stem from deviations between data and the 
OPE approach for the running Rt,v/a{ s q < Mj?) in the vector and axial-vector channels 
(visualized in Fig. 17 of Ref. J7|) and from the corelation between the fitted dimension 
D — 6 and D = 8 operators. They have been found to be larger than the theoretical 
and experimental uncertainties. We repeat this study here in order to estimate the 
corresponding systematic uncertainties for the fitted quantities. The last numbers in 
Eqs. (|29|)-(|3TD, denoted as "OPE" errors, give the deviations found. They are small 
for L\q and dominant for the nonperturbative operators. 

Table fj gives the correlations between the fitted parameters which are found to 
be small. Nevertheless, the interpretation of the parameter errors given in Eqs. 



(pl|) as individual errors must be done with care in the presence of non-vanishing 
correlations. The results can reliably be used when applying the whole expansion ([5]) 
which yields Eqs. (g0|) and fl2j). 

Expressing L\q of Eq. (|29| ) by means of Eq. (|23| ) at the %PT renormalization scale 
/^xPT — 770 MeV, we obtain 

U 10 (M p ) = -(5.13 ± 0.19) x 10~ 3 . (32) 

In deriving the above value the term 2m 2 (25 3 — -£^2,2) in Eq. ( P3"D has been neglected. 
Naive dimensional analysis estimates [|37| give for the low-energy constants B 3 and H 2 .2 
an order of magnitude of 10~ 2 , leading to a contribution which is negligible compared 
to the theoretical error in Eq. (^). The result (|32"D is to be compared with the one- 
loop value of L r 1Q obtainable from 7r — > ewy decays and (r 2 ) n . The error in the latter 



has been decreased since the first determination of L\ Q in Ref. [20]. Using the value 
{r\ = (0.439 ± 0.008) fm 2 ||, one obtains 



L r 9 (M p ) = (6.78 ± 0.15) x 10 _d , (33) 



which updates the value of Refs. pOl |3q|. The L r w constant can be determined from 



the one- loop expression of the axial form factor, Fa, of n — » ez/7 (see Ref. [39| for 
notations): 

Fa = ^^(L 9 + L l0 ) . (34) 



Taking the value F A = 0.0116 ± 0.0016 |3T|, one obtains 

L 9 + L W = (1.36 ±0.19) x 10" 3 , (35) 
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and hence 

L r 10 (M p ) = (-5.42 ± 0.24) x 10~ 3 . (36) 

This is the updated value of L w to be compared with our result (^). Note that the 
quoted errors in Eqs. (|32])-([36]) do not take into account uncertainties from higher or- 
der chiral corrections. A two-loop evaluation for the combination Lg + L\q has been 
attempted in Ref. f3~9fl . However, this analysis is affected by the controversial question 
about the value of the SU (2) constant l 2 . Using for example the values of Ref. |RJ for 
the constants l\ and l 2 on the basis of the two-loop precision tttt scattering phenomenol- 
ogy, we obtain L9 + L10 = (1.07±0.19) x 10~ 3 instead of L 9 + L w = (1.57±0.15) x 10~ 3 
given in Ref. p9|. 



The total, purely nonperturbative contribution to R T y-A found in the fit, taking 
into account the correlations between the operators, amounts to 

R T y- A = 0.061 ± 0.014 , (37) 

compared to the measurement R T y_A = 0.055 ± 0.031. The reduced error of the theo- 
retical fit to data compared to the measurement stems from the additional information 
used in the fit which is obtained from the shape of the spectral functions and the OPE 
constraint. The result fl37|) is in good agreement with the value of R T y — R t ,a = 0.068 
found in Ref. . This is a non-trivial result keeping in mind the logarithmic s de- 
pendence of the dimension D = 6 Wilson coefficients used in this analysis compared 
to the vacuum saturation hypothesis adopted in Ref. @. In addition, in Ref. J7|, vec- 
tor and axial- vector were not combined in a simultaneous fit. The smaller systematic 
error on the nonperturbative parts which is found in this analysis, in particular the 
reduced uncertainty from the explicit dependence of the moments employed, is due to 
the reduced correlation between the fitted D = 6 and D = 8 operators (see Table U). 
The dimension D = 6 contribution to R T y_A corresponding to our fit result Eq. 



amounts to R^y-A = 0-071 ± 0.018, which is significantly than what one obtains from 



the vacuum saturation hypothesis fTJ], R^y^A — 0-97 x 256ir 3 a s (qq) 2 /M® « 0.012. 



In addition to the test of the OPE by varying the (k, I) moments used to fit L 



cff 
10 



and the nonperturbative operators, we perform fits for variable "r masses" so < Mj? 
which provides a direct test of the parameter stability at M 2 . In order to perform such 
a study one has to replace all r masses in Eqs. (0), ( pl|) and (^) by s , while the latter 
must be corrected by the kinematical factor (1 — s/s )(l + 2s/ s )/s . The scale invari- 
ance of the dimension D = 6 operator for variable Sq is approximately conserved when 
keeping the scale parameter [i = M T in Eqs. (H) and (|2l|) unchanged. The dimension 
D = 8 operator is assumed to be scale invariant. Figure ^ shows the fitted observables 
as a function of Sq. The horizontal bands give the results at M 2 within one standard 
deviation. All curves show a convergent behaviour for sq — > M 2 . Any deviation from 
the fitted values for s > Mj? should be covered by the "OPE" errors assigned to the 
results ©-(H). 

Since we use G%PT formulae in this analysis we have investigated the sensitivity 
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1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 



s (GeV 2 ) 

Figure 4: Fit results for Lf^ and the nonperturbative operators as a function of the "r 
mass" Sq. The bands depict the values fl29|)-(|31|) within errors, obtained at M%. 

of the (V — A) t data to a possible constraint on the mass ratio r itself. Clearly a 
combined fit of L\q, r and the nonperturbative operators must fail due to the strong 
correlations of the input variables which reduce the effective degrees of freedom of the 
fit. Thus, as a test, we may use as input for Lf^ and the nonperturbative operators the 
values (PU|)-(PT|) and assume them to be perfectly known, e.g., from a precise second 
measurement. Fig. |5| shows the theoretical prediction of the (most sensitive) IMSR 
moment R^y^\ as a function of r within the errors from the other theoretical sources 
given in Table [l], dominated by the error on a s . Additionally shown as a horizontal 
band are the ALEPH data within experimental errors. We conclude that the current 
experimental precision of the non-strange data does not allow to constrain the light 
quark masses, i.e., the mass ratio r. In the limit of zero u, d quark masses (r — > oo) we 
obtain R^y_ A = 5.11 which is still within the data band of one experimental and the- 
oretical standard deviation. The sensitivity on r when employing the I > moments 
is even worse than with the IMSR. 
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10 15 20 25 30 35 40 45 



50 



r = m /m 



Figure 5: Theoretical prediction of the IMSR moment R^y_ A , using L 
as fixed input value, versus the mass ratio r. The theoretical uncertainty stems mainly 
from the error on a s (M 2 ). The dashed band shows the (V — A) data from hadronic r 
decays (including low energy e + e~ vector cross sections) within experimental errors. 



-eff 



-6.36x10" 



7 Conclusions 

This article deals with a combination of finite energy sum rule techniques and Chiral 
Perturbation Theory (%PT) low-energy expansion in order to exploit recent ALEPH 
data on the non-strange r vector and axial- vector spectral functions with respect to an 
experimental determination of the ;\PT quantity L%q. The theoretical predictions of the 
spectral moments, R^y_ A , of the r hadronic width involve nonperturbative elements 
of the Operator Product Expansion when calculating the contour integral at \s\ = M 2 . 
In the case of inverse spectral moments (I < 0), additional parameters appear 
originating from a second contour integral at the \s\ — 4M^ production threshold which 
subtracts the singularity of the (s/Mj 2 ) -1 inverse moment at s — 0. A constrained fit 
of / < and I > spectral moments adjusts simultaneously the parameter Lf f, defined 
by Eq. (|22|), and nonperturbative power operators of dimension D = 6 and D = 8. We 
obtain Lf f = -(6.36 ± 0.09 ± 0.16) x 10~ 3 , where the first error is of experimental and 
the second of theoretical origin. The present determination of is independent of 
any chiral expansion; in particular, the value obtained here can be directly used in a 
two- loop analysis: it suffices to include higher order corrections in Eq. (^). Within 
the one-loop the above result corresponds to L^ (M P ) = —(5.13 ± 0.19) x 10~ 3 , 
in good agreement with the value L[ (M P ) = —(5.42 ± 0.24) x 10~ 3 extracted from 
the one-loop analysis of ir — > ewy data and (r 2 ) n . This provides a non-trivial test 
of chiral symmetry underlying %PT. The total nonperturbative prediction to i? T y_^ 
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found in the fit, is in agreement with the values of the ALEPH a s (M^) analysis 0. 
The stability of the fit results is investigated in performing various fits for "r masses" 
smaller than M T . Satisfactory convergence is observed. 
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